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Role of Umklapp Processes in Conductivity of Doped Two-Leg Ladders 
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Recent conductivity measurements performed on the hole-doped two-leg ladder material 
Sri4_ x Ca x Cu24 04i reveal an approximately linear power law regime in the c-axis DC resistivity 
as a function of temperature for x = 11. In this work, we employ a bosonic model to argue that 
umklapp processes are responsible for this feature and for the high spectral weight in the optical 
conductivity which occurs beyond the finite frequency Drude-like peak. Including quenched disorder 
in our model allows us to reproduce experimental conductivity and resistivity curves over a wide 
range of energies. We also point out the differences between the effect of umklapp processes in a 
single chain and in the two-leg ladder. 

PACS numbers: 73.50.-h, 72.15.Eb, 72.15.Nj 



I. INTRODUCTION 

Ladder materials have been the focus of much theo- 
retical and experimental work over the last few years jl| . 
These materials, which consist of planes of copper-oxygen 
ladders separated by charge reservoir layers of other 
atoms, are interesting for a number of reasons. In terms 
of structure and behaviour there are similarities to the 
high-T c cuprate materials. For example, hole doped 
two-leg ladders exhibit a spin gap and the doped holes 
form pairs with a d x 2_ y 2 structure. Experimentally, a 
roughly linear temperature dependence of the resistiv- 
ity has been observed in both doped ladders and in 
the high-T c cuprates B. In addition, the copper-oxygen 
"ladders" are weakly coupled to one another, giving rise 
to a quasi one-dimensional system. This is interesting 
because there are analytical techniques for studying one- 
dimensional systems, and, perhaps more importantly, 
one-dimensional electronic systems exhibit Luttingcr Liq- 
uid behaviour characterized by bosonic collective excita- 
tions and spin-charge separation. Ladder systems pro- 
vide an opportunity to study such behaviour, which is 
not normally found in higher dimensional systems. 

Sri4_ x Ca x Cu2404i is an example of an experimentally 
realizable ladder compound This material, in com- 
parison with other ladder materials, is easily fabricated 
and has received considerable attention from experimen- 
talists. The structure of this compound consists of three 
unique layers: a copper-oxygen two-leg ladder layer, a 
Cu02 chain layer, and a layer of Sr and Ca atoms. By 
replacing Sr atoms with Ca atoms, pre-existing holes can 
be moved from the chains to the ladders. For x = 0, the 
ladders are undoped and have a ground state with strong 
antiferromagnetic correlations. 

Conductivity measurements have been performed on 
Sri4_ x Ca x Cu2404i for x — 8 and x = 11 The c- 
axis (parallel to the legs) resistivity, p c , and conductiv- 
ity, (Tic, for x = 11 are of particular interest. For this 
doping the resistivity shows insulating behaviour at low 
temperatures and metallic behaviour above T w 100 K. 
The conductivity has a low-frequency peak located at a 
finite frequency (oj m 50 cm" 1 ). At higher frequencies 



there is considerable spectral weight. 

Conductivity of single one-dimensional chains in the 
presence of either quenched disorder or umklapp pro- 
cesses has been studied in various theoretical papers || ||, 
0, U . The effects of quenched disorder in two- leg ladders 
have also been examined in detail ^ In one such 
paper, Kim [fiof explains the low temperature insulat- 
ing behaviour and the low frequency optical conductiv- 
ity peak of the Sri4_ x Ca x Cu24C>4i two-leg ladder. Also, 
Qin et al. 1 1 1 have recently examined the conductivity 
properties of the doped two-leg ladder using a mean field 
holon-spinon treatment. However, an explanation of the 
metallic behaviour in resistivity above 100 K, specifically 
the monotonic increase in resistivity, for x = 11 has not 
yet been presented. It has been suggested pj that disor- 
der, in the presence of pairing fluctuations, may provide 
an explanation provided spin fluctuations are suppressed 
by the presence of a gap. The purpose of this paper will 
be to examine the roles which quenched disorder and 
umklapp processes play in these regimes. Specifically, 
we demonstrate that disorder cannot be responsible for 
the observed phenomena whereas umklapp processes pro- 
vide a plausible explanation, as they do for the high-T c 
cuprates frL2[ . In addition, we will illustrate the differ- 
ences between umklapp processes in a single chain and 
in a two-leg ladder. 

The plan of the paper is as follows. In Sec. II, we 
present the model Hamiltonian used to describe the two- 
leg ladder with umklapp processes and disorder. This 
Hamiltonian may be derived in the case of weak coupling, 
but becomes a phenomenological model in the strong 
coupling case. In Sec. Ill, we present some scaling re- 



sults which demonstrate that, at high energies, quenched 
disorder gives rise to a decreasing resistivity or a resis- 
tivity with little temperature dependence (depending on 
whether spin fluctuations are present or not) whereas 
umklapp processes can give rise to a resistivity which 
is approximately linear in T. In Sec. [V, we present a 



renormalization group (RG) analysis for umklapp pro- 
cesses and disorder. Using this RG analysis and NMR 
results from Fu jiya ma et al. [jl3| , we argue that the scaling 
results of Sec. [II apply to the Sr3CanCu2404i two-leg 
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ladder. We also discuss how perturbative calculations of 
conductivity, which should be valid at energies well above 
the charge and spin gap energies, are also valid, at least 
qualitatively, until localization effects take over at low en- 
ergies. In Sec. we perform perturbative calculations 
for conductivity using the memory function formalism, 
in analogy with Giamarchi's calculation for the sin- 
gle chain with umklapp processes. Finally, in Sec. [v| , 
we use the perturbative calculations of Sec. [v] along with 
the scaling equations of Sec. Ill to produce conductivity 
curves which resemble the experimental results of Rcf . Q . 



II. MODEL 

A. Clean System 

The Hamiltonian for the clean two-leg Hubbard Ladder 
is given by 



H = -t 



i,s 

f Y (4+1, As +4,s d i+lr 
i.S 

t±_ Y ( c [a,s + 4,s c ^) 

i,s 

:„;.< : -i Y^<-<. ■ (i) 



where the Ci iS 's annihilate fermions in spin state s from 
the i th site on one leg and the di iS 's perform the same 
operation on the other leg. Also, : : represents normal 
ordering with respect to the filled Fermi sea. To investi- 
gate the properties of this model, we begin by making a 
transformation to the bonding/antibonding operators, 



± di 



V2 



(2) 



The resulting free Hamiltonian (U = 0) has a band struc- 
ture consisting of two cosine dispersion relations sepa- 
rated by 2t±, where tj_ is the hopping energy along a 
rung. 

Next we linearize about the Fermi points and take 
a continuum limit as in Ref. The continuum 

Hamiltonian will be composed of chiral field opera- 
tors, ipL/R,b/a,yii where L/R refers to left/right moving 
fermions, b/a refers to the bonding/antibonding band, 
and t / I refers to spin up/down. Specifically, the inter- 
action terms will be of the form e lqx tfjl l ifjgif}s4'~ti with q 
being one of ±{0, 2fc b/a , 4fc b/a , 2(k a ± k b )}, where k a/b is 
the Fermi momentum of the antibonding/bonding band. 
At low energy scales, only the q = terms are important, 
unless one of the other q values is near a reciprocal lattice 
vector (giving rise to umklapp processes). For now we 
will ignore such processes and concentrate on the q = 
case. 



Since we are concerned with the low energy properties 
of the two-leg ladder, a cutoff of 0(l/a), where a is the 
lattice spacing, may be introduced. This allows the well 
known Bosonization transformation | fl5| to be applied to 
the continuum Hamiltonian. Following Lin et al. , we 
perform this procedure and ignore interaction terms with 
irrelevant couplings. The resulting bosonic Hamiltonian 
is 



H — H p + + H p - + H q — h 



where 



and 



dx 



(3) 



(4) 



f dx , — I — I — 

H int — 9c — cos V 47T# p - cos v 2n(j)ba cos v 2nq 



/dx 



g bs COS V 87T0 b(T + g as COS V 87T0 Qcr 



(5) 



A number of comments must be made about this Hamil- 
tonian. First, we have neglected terms proportional to 
the difference between bonding and antibonding Fermi 
velocities in order to eliminate cross terms from the free 
Hamiltonian. This is a good approximation if tj_ is 
smaller than 2t and neither band is very far from half 
filling. Second, in the non-interacting limit, H lrA — > 0, 
u v becomes the Fermi velocity, and the Luttinger pa- 
rameters, {K^}, become unity. Finally, something must 
be said about the boson fields themselves. Within the 
bosonization framework, fermionic densities are repre- 
sented by gradients of boson fields. In the present case, 
charge density in the bonding/antibonding band is given 
by -^d x (j> b/ap while spin density is given by -^d x <j) b / a<7 . 
The remaining fields in Eq. (0) are defined as 



<t>bp ± 4>ap 
4>b<j ± <t)gg 

V2 ' 



and 



7Tj = ~d x 9i 



(G) 



(7) 



Lin ct al. |T6| perform an RG analysis on the con- 
tinuum fermion Hamiltonian for the two-leg ladder. By 
comparing the couplings in Eq. (||) to their counterparts 
in the fermion model, it is readily seen that g c and g b / as 
are relevant in the RG sense. This implies that the two 
spin modes are gapped and one charge mode is gapped. 
Also, in order to pin 9 p - , the two spin modes must first 
be pinned. That is, the magnitude of the charge gap 
must be smaller than the spin gap. Although Eq. (J5|) 
may be derived from the Hubbard for small U , numeri- 
cal work [jl7| shows that there is no phase transition as 
one moves from weak to strong coupling. Therefore, we 
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take this Hamiltonian as a phenomenological model for 
the case of strong coupling as well. 

The free part of the Hamiltonian in Eq. (S) is written 
in terms of the (j> a ± fields to eliminate cross terms which 
would appear if it had been written in terms of the 4> b i aa 
fields. The interaction part of the Hamiltonian is written 
in terms of the b/aa fields in order to elucidate the gap 
properties of the system. However, it is convenient to 
express the entire Hamiltonian in terms of the cr* fields 
when describing scaling and renormalization properties. 
The interaction part then becomes 

Hint = ^\ J dx cos VAnOp- cos ^JTn4> a + 

92 



dx cos v An9 p - cos v 47T</) cr - 

+ — £ / dx COS V47T0CT+ COS \f^K<j> a - 

a z J 

+ —5- / dx sin V47n/> (T + sin VAtt^o— 1 (8) 
a 1 J 

where g\ and g 2 are equal to g c , but will flow differently 
under RG, and g± = (g as ± 5fc s )/2. 

B. Scattering Mechanisms 

1. Umklapp Processes 

In the previous section, interaction terms of the form 
e^t/j^ip^sip-f were ignored for non-zero q values. How- 
ever, if any of the possible q values are close to a recip- 
rocal lattice vector, then such terms will become impor- 
tant. With the restrictions on t± and filling (t± small 
compared to t, and both bands not far from half- filling), 
the only q values which can satisfy this condition are 
±{4fcf,/ a , 2(fc a + k b )}. In these cases the net momentum 
transfers will be 2tt /a — q. When Bosonized, the umklapp 
terms become 

/dx r 
— cos (V47r0 p+ + S b+a x) g u i cos V4tt9 p - 



g u2 cos V 4:TT(f) a - - g u3 cos V AnO, 



dx 
as 
dx 

g„r, I — cos 



47T I 



47T I 



S a x 
S b x 



where 



and 



2tt 

h+a = 2(fc a + k b ), 

a 



X 2n AT 

hi a = 4/e 6/a . 

a 



(9) 



(10) 



(11) 



Eq. (Q) reveals an important difference between umk- 
lapp processes on a single chain and on a two-leg ladder. 



For the chain, spin was completely decoupled from the 
charge carrying mode jjj. For the ladder, umklapp pro- 
cesses couple the charge mode, p + , to spin. The conse- 
quences of this coupling will be discussed later. 



2. Disorder 

For one dimensional conductors, disorder has typically 
been modeled by a random on-site potential p, 0, |l8fl, 



(12) 



where p is the leg label, i is the site label and s is the 
spin state label. Also, the disorder strength, e^ )P , is taken 
from a Gaussian distribution of width, D, which should 
be proportional to vf/t, where vf is the Fermi velocity 
and r is the mean free path. 

After switching to bonding/antibonding operators, 
taking the continuum limit and Bosonizing, the disorder 
Hamiltonian becomes 



Ha 



dx 



— t, 2 e v yv p + 



dx 

T7« 



' ' cos \fir(<j>, 
p~ ' cos \/tt{4>, 



+ 



■) 



e 
H.c 



+e p - ) cos (0 ct+ + e a _ )) 

i ^p-+ e p-)cosV^(0 ff 
iv^(0 P --e p -) cos ^ ( ^_ _ 



IX' 



>) 



(13) 



where the £'s are taken from independent Gaussian dis- 
tributions. Specifically, the disorder fields satisfy 



((^{x)^(x'))) = aDiSijSix - x'), 

«££» = 



(14) 



where ((...)) denotes disorder averaging. Terms which 
produce scatterings from a Fermi point to itself have been 
ignored, as they are unimportant for conductivity cal- 
culations and may be eliminated by simple transforma- 
tions M. 

Eq. (|13j) is useful for performing perturbative calcu- 
lations above the spin and charge gaps. However, be- 
low these energies, simplified forms for the Hamiltonian 
may be derived by integrating over quantum fluctua- 
tions which arise due to the pinning of a field and the 
rapid oscillations this induces in the corresponding con- 
jugate field. This procedure is explained in Appendix C 
of Ref. §. 

Some earlier experimental evidence [ [l9| points to a 
spin gap of about 300 K for Sri4_ x Ca x Cu2404i, with 
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x = 11.5. One possible explanation of the nearly linear 
power law dependence of p c on T JTcf] , as mentioned in 
the introduction, involves quenched disorder in the pres- 
ence of pairing fluctuations without spin fluctuations (i.e. 
Op- is free to fluctuate, but <fi a + and <p a - are pinned). If 
such a regime exists, the effective Hamiltonian may be 
derived by integrating over the fluctuations in and 
9„- . The result is 



H 



dx 



dx 

-i\A r (^„+ -<t>„- ) 



H p + + H p - +g c I — cos V4ir6 p - 



6e 



9-> +H.C. , 



(15) 



where terms which would contribute to physical proper- 
ties by an amount of 0(D 2 ) have been dropped. 



III. CONDUCTIVITY AND SCALING 

The real part of the conductivity is given by 
-1 



-9f[n(w,T)] 



(16) 



where II is the current-current correlation function which 
may be calculated using the Kubo formula. In one dimen- 
sion we have 



n(w) 



dx / dte^{\p\x,t),3$,Gj\), (17) 
Jo 



where j is the current operator, which can be readily 
obtained from the charge density operator, p, using the 
continuity equation, d t p — d x j = 0. For the present sys- 
tem, p is given by 



P 



-J) 



x<Pp+, 



yielding 



--dt<t>p 



Inserting this into Eq. (17) gives 



n(w,T) = — G Iet (cj,T) 

IT 



(18) 



(19) 



(20) 



where 



G rct = / dx dt e iwt ( [0 p+ (x, t), <P p+ (0, 0)] ) . (21) 

At high energies, the retarded Green's function, G re t, 
may be calculated using perturbation theory. However, 
if low energy properties are of interest, then another ap- 
proach must be taken. In this work, RG and scaling 



techniques [ f20[ are used to obtain low energy informa- 
tion from perturbative calculations of higher energy prop- 
erties. Suppose energy is scaled in the following way: 
E -> bE (b = e dl and E = uo or T), then the Green's 
function transforms as 



G let (E, {g}) = b 2 Z p+ G let (bE, {g(dl)}) 



(22) 



where Z p + is due to "wave function renormalization" , {g} 
represents the set of coupling constants, and {g(dl)} rep- 
resents the set of coupling constants after scaling. Letting 
Z p + = b 1 and iterating Eq. (E2) gives 



G Iet (E,{g}) = exp 



dl (2-7(0) 



x Grot (e r E,{g(l*)}). 



(23) 



If we wish to calculate the Green's function on the LHS 
of Eq. (p3]), we may choose 1* such that e l " = (E /E) 
where E a is an energy comparable to the bandwidth (i.e. 
perturbation theory should be valid there). 

Eq. ( ^3|) is also useful when E is large and the Born 
approximation is valid. In this case, it leads directly to 
a scaling relation for 1/r. It is 



(24) 



-(bE,{g(dl)}) r(E,{g})' 



where contributions of 0(g 2 ) have been ignored. Since 
resistivity is proportional to 1/r, Eq. ( p4| ) may be used 
to calculate its temperature or frequency dependence in 
regimes where perturbation theory is valid (i.e. well 
above spin and charge gaps). In order to do this, the 
scaling behaviour of couplings which lead to scattering in 
the p + channel must be known. This information is eas- 
ily determined in the Boson formalism (see Ref . |H| ) . As 
an example, the disorder strengths, D a / b , in the Hamil- 
tonian of Eq. (§||) scale as D' a/b = D a/b b 3 - {K e+~ K p- )/2 . 
This yields a high temperature resistivity of the form 



p c (x a x D a T 2 



a 3 D b T— 



(25) 



where the a^s are constants. This form for resistivity is 
valid above the p~ gap if spin fluctuations are sufficiently 
suppressed. In the next section, we argue that the Lut- 
tinger parameters cannot deviate significantly from unity 
without moving the spin and/or charge gaps to extremely 
high or low energies. Therefore, Eq. (25) implies that p c 



will have a monotonically decreasing temperature depen- 
dence in the presence of pairing fluctuations and disorder. 
This is inconsistent with the linear power law regime ob- 
served in Sr 3 CanCu2404i. 

Although some earlier experimental evidence points to 
a 300 K spin gap for the strontium cuprate ladder, NMR 
results from Fujiyama et al.[]l3| indicate the existence of 
spin fluctuations at temperatures as low as 60 K - about 
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the same temperature at which charge localization oc- 
curs. The Hamiltonian consisting of Eq. (|J) and (13), 
which allows for spin fluctuations, may, therefore, be a 
more appropriate model in the power law regime. Per- 
forming the same scaling analysis which led to Eq. ( p5| ) 
on this model yields a high temperature resistivity of the 
form 



a 3 D 2 (£) 
a 5 D 3 (f ) 



2 -Kv+ A V+^++ K »-) 



2-i( K ,+K _+K 



/a 8 \ 2 -H K P++ K P-+ K *-+-i 



(26) 



For Luttinger parameters close to unity, this resistivity 
will have little temperature dependence. 

Turning to umklapp processes, if the system is not far 
from a special filling (i.e. one of the 5's in Eq. (^) is 
small), then, at high temperatures, Eq. (|24|) can be used 
to determine p c as before. For the case of 8b+ a small, it 



2 ( a 2 \ 3 ~ 2 ( K p++TT 



Pc oc aig ul y 



2 /04 

m 



3-2( V + TTT 



(27) 



If the Luttinger parameters are close to unity, then this 
resistivity is approximately linear in T . For the case of 
5b/ a small, similar results hold. 

At this point, an important difference between the ef- 
fect of umklapp processes on the conductivity of single 
chains and two-leg ladders becomes apparent. In the case 
of a single chain, if 8 <§; T but not zero (8 = gives rise to 
a charge gap) and K p w 1, then umklapp processes lead 
to a linear T dependence for resistivity regardless of the 
behaviour of the spin mode (which is completely decou- 
pled from the charge mode) |7j. However, in the case of 
the two-leg ladder, charge and spin gaps can alter this be- 
haviour. For example, if a specific ladder material has a 
large spin gap, then a regime in which fluctuations of the 
p~ mode are allowed but spin fluctuations are suppressed 
might exist. For 8b+ a small, an effective Hamiltonian for 
this regime may be derived from Eq. (^J) by integrating 
over fluctuations in the Q a - field. The result would be 
terms in the DC resistivity which are proportional to T 5 , 
T and I/T. Below the p~ gap, all fillings which give rise 
to significant umklapp scattering would give T 5 and 1/T 
contributions |2l|. 



IV. RENORMALIZATION GROUP ANALYSIS 

In the previous section it was demonstrated that, at 
high energies, umklapp processes can contribute a linear 
temperature dependence to the DC resistivity of two- 
leg ladder compounds. In this section and the next, we 
discuss what happens at lower energy and attempt to 
explain the conductivity data of Osafune et al. 0] Based 
on the discussion of the previous sections, we take our 
model to consist of the clean Hamiltonian of Eq. (Bj) , the 
8f, +a term of the umklapp Hamiltonian of Eq. <M) (we 
drop the b + a subscript from here on) , and the disorder 
Hamiltonian, 



H, 



dis 



dx 
a z 



— i v4?r0 



(28) 



used in Rcf.[|10|. Note that Hdis can be derived from 
Eq. (13) by integrating over fluctuations in the and 
the p~ modes which exist when g^/as and 9c are large. 
Since these couplings are strongly relevant, Eq. (|28| ) is the 
correct disorder Hamiltonian at low energies. Here, we 
assume that the disorder couplings appearing in Eq. psj) 
are sufficiently weak that H^is has little effect on trans- 
port properties except at these low energies. 

RG equations may be derived for our model by scal- 
ing space-time variables while adjusting coupling con- 
stants to maintain the real space cutoff, a. Following 
standard methods ^0|, this procedure has been carried 
out to second order in the couplings, with the results 
shown in Appendix |X| These equations have a number 
of important implications. In weak coupling, the and 
the p~ modes are gapped due to the relevancy of gi, g 2 , 
g as and gb s - Numerical work fll7| has shown that these 
gaps persist in the strong coupling case. Therefore, the 
parameters of our phenomenological Hamiltonian must 
have values such that the couplings which generate the 
gaps grow large at the appropriate energy scale. From 
the RG equations it is also clear that, when the inter- 
action couplings are small, the energy scale of the gaps 
is determined by the Luttinger parameters. For the case 
of a one-dimensional chain, there is only one Luttinger 
parameter, K a , relating to spin. This parameter is re- 
stricted to values near unity |22| . We will assume that it 
is reasonable to take K a ± near unity as well. In fact, if ei- 
ther of K a ± or K„- deviates significantly from unity, then 
one or more gaps will open at very high energy (near Eq) 
or at very low energy. This would be inconsistent with 
a gap magnitude of approximately 60 K. Taken together 
with the fact that K p + < 1 for repulsive density inter- 
actions, these Luttinger parameter restrictions indicate 
that umklapp scattering gives rise to a nearly linear T 
dependence for resistivity at high energies (see Eq. (p7|)). 
In addition, numerical integration of the RG equations 
shows that the Luttinger parameters do not undergo sig- 
nificant renormalization until the gap energy is reached. 
Therefore, this linear T behaviour should persist until 
just above the gap energies, well below Eq. 



6 



Given that scaling arguments strongly suggest umk- 
lapp processes as the mechanism responsible for the 
power law resistivity feature of Sr 3 CanCu2404i, it would 
be useful to perform conductivity calculations which are 
valid at all energy scales. To accomplish this, pertur- 
bative calculations along with the RG equations and 
Eq. ( p2] ) may be used. A number of simplifications can be 
made at this point. First, from numerical integration, we 
find that the Luttinger parameters and the velocities flow 
slowly at high energies (their flow rate is proportional to 
0(g 2 )). At low energies, renormalization due to umklapp 
processes is attenuated by a rapid growth of S. There- 
fore, we ignore renormalization of Luttinger parameters 
and velocities due to umklapp processes. Within this ap- 
proximation, the velocities will all flow to the same value 
until localization takes over at low energy. This implies 
that, above the localization energy, there is no qualitative 
difference between conductivities calculated with differ- 
ent velocity parameters. Taking all velocities as equal 
also allows perturbative conductivity calculations to be 
performed in a straightforward manner using standard 
results f 
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Finally, all gap generating couplings 
will flow slowly under renormalization until they become 
0(1), at which point they diverge rapidly. The Luttinger 
parameters also flow slowly until this point. Given that 
the Sri4_ x Ca x Cu2404i ladder has spin fluctuations at 
energies as low as the charge localization energy, we ig- 
nore renormalization of the Luttinger parameters. This 
is qualitatively valid since these parameters only flow in 
a significant way after wavefunction renormalization be- 
comes the dominant factor in scaling (see Eq. (B2n). 



V. PERTURBATIVE CALCULATION OF THE 
MEMORY FUNCTION 



In order to make use of Eq. ( p2| ) , perturbative calcula- 
tions of p c and a c are required. Following Giamarchi |Q, 
the memory function formalism |^5[ may be implemented 
for this purpose. We proceed by writing the conductivity 



a = \Q- 



M(uj) 



(29) 



where Q is the response of the free electron system and 
M is the memory function given by 

MM = <^oH-(fi (30) 



ujQ 



with 



oc poo 



(F; F) Q = -i / dx dt 4 ([F(x, t), F(0)]) 

-oo JO 



(31) 



and F = [j, H}. From Eq. (§), @ and (§|) we find 



F = 



i cos 



- g U 2 cos \/47r</> -- ~ 9u3 cos V^r6 a - j 
+ r(z)e ix/i¥0 ' + -£(a;)e- i73 ^ 



4tt6> 



(32) 



where only the Sb+ a term from the umklapp Hamiltonian 
has been used and all velocities have been taken as equal. 
Using the results of Schulz |2^j] , and Schulz and Bourbon- 
nais P4| , (F; F) (ui) may be calculated in a straightfor- 
ward manner. The memory function of Eq. (|3^) then 
becomes M — M um + Afdis, with 



and 



ttK p + 1 „ , /27raT\ 2(Kx_1) 



(33) 



Ma 



AD sin(7rA" 



/27raT\ 



(34) 



where 



f-isAl-K p 



x B 



B >0 



= B 



K 



U) ± ux 



and 



4ttT 



K„- 



X 



1 

X = 2 
X = 3 



(35) 



(36) 



Here B(x,y) is the Beta function and Q has been taken 
as 4uK p + /n, twice the value for a single band model p6[ . 
Also, the p integral in Afdis arises because the disorder 
field, £, is Fourier transformed before disorder averaging 
is performed. 

In this work, the zero frequency and zero temperature 
limits of the memory function will be of interest. For 
T^O.w^Owe have 



M, 



ttK, 



[(u6) 2 -oj 2 X -(uS) 2(Kx 



2(«x"l) 



(37) 
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and 



M dis -> 
4DK n + 1 



^ sin(7rif p+ )r 2 (l- J ft: p+ 



(up) — of* 



■K.4--1 



(up) 



I \ 2 



(38) 



For the case u — > 0, T ^ we have 

3 



M„ 



— 12 - 2^ 3n x r (1 - K X ) Sm(TTK x ) 

x=i 

f2-KaT\ 2( - Kx ' 1} f 1 \ T(z x )r(z x ) 



and 
M dii 



\ u J \a-kT J r(i - %)r(i - z x ) 

x [*(1 - z x ) + *(1 - z x ) - *(z x ) - (39) 



i4Di ^ir 2 (i-if p+ )sin(^ p+ ) ^ — 



(7^ 



r (y P ) r (y P ) 



/ 27rar \ 2( ^ + - 1) 

\~^~) J^ r " e T(l-y p )T(l-y p ) 
x [tf (1 - y p ) + *(1 - y p ) - - $(y p )] (40) 

where *(a;) = r'(a;)/r(x), 

K v . u5 



and 



2/p 



4ttT 



K„+ . lit) 
— r- +1 — — 



4iraT 



(41) 



(42) 



VI. RESULTS FROM SCALING 

Using the memory functions fro m the previous section 
along with Eq. ( |l6| ) , (|||) and (A.1), conductivity and re- 
sistivity curves may be produced for a given set of cou- 
plings. To proceed, it is useful to rewrite Eq. (^9| ) using 
Eq. (16) and ( p2| ) to obtain a scaling relation for resistiv- 
ity/conductivity in terms of the memory function. For 
resistivity we have 



P(T, {<?}) = 



AuK p+ (l*) 



Jo dl ^ l )M(e l 'T, {g(l*)}) 



(43) 



and for conductivity, 



$ta(u,{g}) 



e l *uK p +(l*) 



°J*uJ + M(e l *uj,{g(l*)}) 



In Sec. IV, it was argued that the Luttinger parameters 
and velocities do not renormalizc significantly until the 
charge localization energy is reached. Therefore, in us- 
ing Eq. (|43| ) and (44), renormalization of these quantities 
will be ignored. Within this approximation, the param- 
eters which affect the resistivity / conductivity curves are 
D, g U i,u,6, K p +, K p - and K a - . Given any set of these 
parameters, we can now generate such curves. To set 
units, we make the transformation 



T- 


-> k B T, 


LO - 


■+ hw, 


9i ~ 


-> g^a, 


u - 


-> uta, 


D - 


■+ Dt 2 a 2 , 


p - 


1 

—Pi 


a - 


+ e 2 a 







(45) 



(44) 



where e is the electron charge. 

By trial and error, it is simple to choose parameters 
which generate curves that resemble those measured by 
Osafune et al. || We first recall that the Luttinger pa- 
rameters should not be far from one, Eq must be well 
above the highest temperature or wavenumber of inter- 
est, K p + must be less than one and 6 should be small. 
It is also of interest to note that, quantitatively, the re- 
sistivity and conductivity curves are most sensitive to 
changes in D and K p + . For both resistivity and conduc- 
tivity, we choose Eq = 0.25 eV, ut = 1.1 cV, K p + = 0.9, 
K p - = 1.15, K a - = 0.85, g uX t = 0.4 eV, g u2 t = eV, 
g u3 t = 0.5 eV, S = 0.022 and D = 0.00034 eV 2 . The 
resulting resistivity curve, shown in Fig. [II demonstrates 
insulating behavior beginning below 60 K and a mono- 
tonically increasing power law regime above that tem- 
perature. The localization temperature is mainly deter- 
mined by K p + , the exact nature of the power law above 
60 K is determined by the Luttinger parameters, the sep- 
aration between the two regimes is determined by 8 and 
the relative heights of the various parts of the curve are 
determined by {g U i} and D. The velocity, u, plays a role 
in determining the overall scale of the curve, as well as 
how fast the resistivity diverges at low T . Our conduc- 
tivity curve, shown in Fig. ^|, demonstrates a finite fre- 
quency, Drude-like peak with additional spectral weight 
at higher wavenumbers. The position of the peak is deter- 
mined mainly by K p + , the peak height and attenuation 
are determined by u and D, the relative height of the 
spectral weight beyond the peak is determined by {g U i}i 
the nature of the power law at high wavenumbers is deter- 
mined by the Luttinger parameters, and the wavenumber 
at which absorption due to umklapp processes begins is 
determined by 5. 

Although we have not done a detailed numerical fit 
to the data, our resistivity and conductivity curves show 
good qualitative agreement with the measurements of Os- 
afune et al. ||. However, it it also clear that one can 
achieve a better quantitative fit to the data by choosing 
different values of D and A^ p + to independently fit the re- 
sistivity and conductivity curves. This quantitative dis- 
crepancy is likely due to the fact that the renormalization 



8 



a 4 

E 



100 200 
T(K) 



300 



ladder compound Sri4_ x Ca x Cu2404i can be explained 
by a combination of disorder and umklapp effects, ft 
has been shown that quenched disorder, whether in the 
presence of pairing fluctuations alone or in the presence 
of pairing and spin fluctuations, does not contribute a 
monotonically increasing temperature dependence to the 
resistivity at high T. However, if spin and pairing fluc- 
tuations persist down to the charge localization temper- 
ature, then umklapp process can contribute an approx- 
imately linear temperature dependence if the system is 
near an appropriate filling. Experimental evidence from 
Fujiyama et al. suggests that such fluctuations do 
exist. 



FIG. 1: p c (T) vs. T. 



APPENDIX A: RG EQUATIONS 
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FIG. 2: oi c (u) vs. wavenumber. 



of the Luttingcr parameters and velocities was ignored. 



We have derived RG equations for the model described 
in Sec. IV using standard techniques. To second order in 
the interaction couplings they are 



^=91 2 



92 2 



dgi 

dl 

dg 2 
dl 

dgbs 
dl 

dg as 
dl 

dg ui 
dl 

dg U 2 

dl 
dg U 3 

dl 



K a - , 



g bs (2-K a+ -K a -), 

g as {2-K a+ ~K a -), 
I 



g u i 2 - K + - 



K„ 



gui (2-K p+ -K a -] 



g u3 I 2 - Kpi 



VII. CONCLUSIONS 



In this paper, we have used RG and scaling techniques 
to argue that the conductivity properties of the two-leg 



dl 



dl 
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du a - 
dl 



dK p+ 
dl 

dKp_ 
dl 

dK a+ 

dl 
dK a - 

dl 



8u 



- 1 + 



{92 + <?+ + 9- + 2J (6a)g 2 u2 ) K a - + 



2 Jo(Sa)g 2 u3 \ 
K a - J 



2J 2 {5a)g 2 u2 K a - + 



2MSa)9l 3 
K„- 



(2ir) 2 (g 2 ul +g 2 u2 + g 2 u3 )K 2 



l + \ M 5a)+ 1-^ )j 2 (5a) 



SttDK 2 + 

U p +Uq 



+ 




g 2 + g 2 + 2J {6a)g 2 ul 



(g 2 +g 2 + +g 2 _)K 2 



M Sa )9li 



(g 2 2 + 9 2 + + 9 2 - + 2Jo(Sa)g 2 u2 ) K 2 a _ - 2J Q (5a)g 2 u3 j 



2J 2 (Sa)g 2 u2 K 2 a . - 2J 2 (6a)g 2 u3 



(Al) 



where J / 2 are the zeroth/second order Bessel functions, 
the b + a subscript has been dropped from 5b+ a , and u$ 
was introduced to set the units of the time cutoff cor- 
rectly. It is immediately clear that the velocities and 
Luttinger parameters initially flow slowly since their rate 
of change is 0(g 2 ). Also, 5 grows exponentially so that 
the Bessel functions which appear in the equations for 
the velocities and Luttinger parameters go to zero at low 
energies. The contribution of umklapp processes to the 
renormalization of velocities and Luttinger parameters 
is neglected in this work by dropping all terms involv- 
ing Bessel functions. Within this approximation, the 
velocities other than u p + will all flow to uq while u p + 
will only flow in a significant way when localization takes 
over at low energies. At this point, wavefunction renor- 
malization becomes the dominant effect so we may ig- 
nore renormalization of the velocities as well. Finally, 
the feedback between the Luttinger parameters and the 



couplings which renormalize it is such that the Luttinger 
parameters flow slowly until couplings become 0(1). At 
this point, both couplings and Luttinger parameters di- 
verge rapidly. However, if the gaps in the system oc- 
cur below the charge localization energy, then this won't 
happen until wavefunction renormalization becomes the 
dominant effect. Under these circumstances, neglecting 
Luttinger parameter renormalization should not intro- 
duce any qualitative errors in calculations of conductivity 
(although quantitatively there will likely be some effect). 
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